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We show that in anisotropic or layered superconductors impurities induce a van der Waals 
attraction between flux lines. This attraction together with the disorder induced repulsion may 
change the low B - low T phase diagram significantly from that of the pure thermal case considered 
recently by Blatter and Geshkenbein [Phys. Rev. Lett. 77, 4958 (1996)]. 
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Conventional type-II superconductors show in addi- 
tion to the flux repulsing Meissner state a second su- 
perconducting (Abrikosov) phase in which the magnetic 
induction B enters the material in the form of quantized 
flux lines (FLs) which form a triangular lattice. Each FL 
carries a unit flux quantum $0 = hc/2e. The Abrikosov 
lattice is characterized by a non-zero shear modulus cge, 
which vanishes at the upper and lower critical field, 
and TJci , where continuous transitions to the normal and 
the Meissner state, respectively, occur. In his mean-field 
solution Abrikosov treats FLs as stiff rods. Close to the 
lower critical field iJci , their interaction becomes expo- 
nentially weak and hence the FL density l^'^ = -B/$o 
vanishes as |ln/i|~^ where h = {H — Hc^)/Hci denotes 
the reduced field strength 

Thermal fluctuations roughen the FLs resulting in a 
possible melting of the Abrikosov lattice close to iJci and 
Hc2, respectively, because of the softening of cge- This 
applies in particular to high-Tc materials with their ele- 
vated transition temperatures and their pronounced layer 
structures Q| . At present, it is not clear whether the tran- 
sition to the normal phase at high field happens in these 
materials via one or two transitions. However, melting of 
the FL lattice has clearly been observed experimentally 

i- 

At low fields a first order melting transition to a liq- 
uid phase and a change in the critical behavior of B has 
been predicted some time ago by Nelson Q. Quantita- 
tively the influence of thermal fluctuations is described 
W a thermal length scale Lt = ^1/{1&tt'^T) « 2cmi^/T 
Lt has a simple physical meaning: an isolated flux 
line of length Lt shows a thermal mean square displace- 
ment of the order of the London penetration length A. 
Besides a shift of TJci, large scale thermal fluctuations 

lead close to Hc^ to an entropic repulsion (^X^/LtIY 
between FLs which dominates over the bare interaction 
for small h and hence ~ /i (Here we measure all 
FL interactions in units of £0 = (*i'o/47rA)^ = LtT/X^). 

More recently. Blatter and Geshkenbein ||] found that 
in anisotropic or layered superconductors short scale fluc- 
tuations give rise also to an attractive van der Waals 
(vdW) interaction 0. For FLs separated by a dis- 
tance R the strength of this interaction is of the order 
of \^/LT{d + eR)R'^. = m/M 1 denotes the 



anisotropy of the material with m and M the effective 
masses parallel and perpendicular to the ab plane, and d 
the interlayer spacing. A and X/e are then the screening 
lengths parallel and perpendicular to the layers, respec- 
tively. 

The competition among the bare, the entropic and the 
vdW interactions leads to an interesting phase diagram 
at low B values. In particular. Blatter and Geshken- 
bein 1^ find at low T a first order transition between the 
Meissner and the Abrikosov phase. 

So far fluctuation effects have been discussed for a 
clean superconductor. It is well known, however, that 
in type-II superconductors FLs have to be pinned in or- 
der to prevent dissipation from their motion under the 
influence of an external current. Therefore, besides the 
thermal fluctuations one has to take into account the 
effect of disorder. Randomly distributed pinning cen- 
ters lead indeed to a destruction of the Abrikosov lattice 
but as has been recently shown, for not too strong 
disorder FLs may form a (Bragg-) glass phase which is 
characterized by quasi long-range order of the FL lattice 
For low B this phase undergoes a melting transi- 
tion to a pinned liquid state. Inside this phase, disorder 
induced effects are expected to dominate over those of 
thermal fluctuations for sufhciently low T. It is therefore 
the aim of the present paper to consider the influence of 
the disorder fluctuation induced forces between the FLs, 
and to study the low B phase diagram. In particular 
we find, that the latter deviates substantially from that 
found in Ref. |^ for pure systems. The rest of the paper 
is organized as follows. We first reconsider the disorder 
mediated steric repulsion between the FLs and then de- 
rive the disorder induced vdW interaction. Finally we 
discuss the phase diagram at low B and T. 

Let us start with the steric repulsion, which results 
from the long wave length fluctuations of FLs. Using a 
simple scaling argument, it was found in Ref. [Tc| ], that 
the disorder dominated steric repulsion between FLs is 
of the order 



(1) 



where the roughness exponent C of a single FL in a ran- 
dom potential is about 5/8 in I? ~ 3 dimensions [ pT| . 
Here ^ is the correlation length in the superconducting 
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planes and Ldis denotes a disorder-related length scale 
with a similar meaning as Lt- At low temperatures, 
Ldis ~ LtTk^/'^~'^ /Td\s, where k = A/^, denotes the 
Ginzburg-Landau parameter. Tdis is defined in (||), and 
it is identical with T^^° of ref. The steric interac- 
tion results from the fact, that the FL is confined to a 
cylindrical cage, formed by its six nearest neighbors. It 
therefore cannot gain the full energy decrease a free FL 
can obtain from the disorder. 

This argument has been later put into question [ p^ , 
since the FL might overcome the averaged repulsive po- 
tential of the cage and hence leave it to fol- 
low its optimal path. Naively, this should be the case if 
(A/O^-i < A2/CL-2;2-2/c or D > (2/CP)) - 1- 
Such an argument is indeed correct for thermal fluctua- 
tions, where C = 1/2 in all dimensions: for Z? > 3 the 
steric repulsion becomes ineffective since the entropy gain 
from leaving the cage outweighs the repulsion fl^ . 

In the case of disorder, the situation is however differ- 
ent. At low T a FL which leaves its cage can only increase 
its energy gain if it follows the optimal path. Let us ne- 
glect for the moment the repulsive interaction between 
the lines. Since the optimal path for a line with a given 
initial point has a river delta like shape, two lines sepa- 
rated initially by a distance nl will merge after a tran- 
sient region of length {niy/'^ <^ L. Let us now switch 
on again the repulsion. Since the lines follow the same 
path, repulsion leads to an energy increase of the order 
L. This has to be compared with the free energy gain by 
following the optimal path instead of the path inside the 
cage, which is of the order i^^^^ and hence smaller than 
the repulsion since C < 1. Thus, a line leaving its cage 
would increase its energy instead of decreasing it unlike 
the thermal case, where lines can avoid each other and 
still gain entropy. Consequently, the simple argument of 
Ref. applies which results in (^. 

The above argument can indeed be derived in a some- 
what more formal way by using a renormalization group 
(RG) analysis of interacting FLs ||l^. In the large wave- 
length limit, the FLs of length L are described by the 
Hamiltonian 



n= / dzV{^R2+£pi„(R„ 



z) + Y,vo5x{R,,)}, (2) 



where Ri = dR.i{z)/dz and Ry = Ri — Rj. (R.i{z),z) 
denotes the position vector of the i-th FL, ei its stiffness 
constant which is of the order eq in the lon g waveleng th 
limit. The random pinning potential fulfills £pin(R, z) = 
and § 



£pin(R,z)epi„(O,0) = (T^Jeon '5(z)fc(R/0 



(3) 



with fc(x) = 1 for X <C 1 and fc(x) w {l/x'^)lnx for 
X ^ 1, respectively. T^j^ is proportional to the impurity 
concentration. (Sa(R) represents a (5-function, smeared 
out over a scale A and vq w AttX^Sq. Since we are inter- 
ested in the diluted limit I :s> X, the RG flow of A and 



V can be obtained by considering only two interacting 
lines as in Ref. (wj. Close to Z? = 2 dimensions, we find 
(dv/dl) — 2(1/C — l)v and hence v remains a relevant 
perturbation even above two dimensions, in contrast to 
the results obtained in Ref. We, therefore, conclude 
that the result (1) is valid for all D. It is worth remarking 
here, that in contrast to the present case, a single FL in a 
cylindrical cage does indeed show a depinning transition 
from this potential in > 2 dimension |]l5|. 

Next we consider the attractive van der Waals interac- 
tion. We start with the case of extreme anisotropy £ = 0, 
in which the Josephson coupling between the layers can 
be neglected. The repulsive interaction between two pan- 
cake vortices 1 and 2 separated by a distance R in the 
same layer is then larger than their attractive interaction 
in different layers by a factor X/d (3> 1). In consider- 
ing the interaction of different FLs consisting of pancake 
vortices we therefore restrict ourselves to pancakes in the 
same layer. Since their interaction is 2c?ln(^/i?), the in- 
teraction of two pancake dipoles resulting from the dis- 
placements Ui , U2 of pancakes 1 and 2 in the same layer 
is Ui2 — — (2(i/i?2)[2(ui • ri)(u2 • n) — Ui • U2], where n is 
the unit vector in the direction connecting the two FLs. 
A displacement Ui will indeed induce a displacement U2 
because of the force fi2 = —dUi2/d\i2. The actual re- 
sponse of the pancake 2 is limited by the elastic force 
f22 ~ — (du2/A^) \a{Xn/d) resulting from the other pan- 
cakes of FL 2. Here we have used the stiffness constant in 
the decoupled limit §, £,(fc) = {eo/2X^k'^) ln(l + X^P) 
for fc « Tr/d. With fi2 + £22 = and averaging over dif- 
ferent configurations of Ui , the vdW-interaction per unit 
length is 
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Inserting the result for the short wavelength thermal fiuc- 
tuations (u^) = {2T/eod) [X"^ / \n{nX/d)], we obtain the 
vdW attraction given in Ref. |^ , apart from a numerical 
factor. 

In the case of an impure superconductor (u^) is given 
by the mean square displacement of a pancake in a po- 
tential which is a superposition of a parabolic elastic 
part £od(u/27rA)2 In {Xir/d), resulting from the dispersive 
elastic constant for k « 7r/d, and a random potential 

(iz£pin(u, z). Rewriting (H) for L w 0? in dimensionless 
quantities r — R/^ and t — z/d, it is easy to see, that 
the ground state displacement Uq/^ is only a functional 
of (A(7rA/d))^/2£pin(r, t). £pin(r,i) is the dimensionless 
random potential of mean zero and variance unity and 
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Thus at r < Tdis, (u ) fvul = C^/(A(7rA/d)). Perturba- 
tion theory applies for A <C 1 and gives /(A) « A** with 
rj — 1. For larger A no exact result is known. Imry-Ma 
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arguments and a variational treatment give 77 = 1/2 with 
a logarithmic correction fl^ , whereas we find 77 w 4/9 
from simulations. For BSCCO with e = 1/300, A w 
2OOOA, C ~ 20A, d = I5A, and Tdis = 45^ |l, Aq « 91 
and A « 261, i.e. we are in the non-perturbative regime. 
In Ref. jl^ a much larger value of Tdis = 21 OX was re- 
ported which results Ao ~ 9.261 x 10^ and A « 2.6 x IQ-*. 
The final form for the disorder induced vdW interaction 
in the decoupled limit is therefore 



V 



(dis) 



vdW 



In TrA/d 



Ao 



TlX/d 

\n{T:\/d) 



(6) 



From (^)- (H), we conclude that, in the decoupled limit, 
the vdW interaction will be dominated by disorder fiuc- 
tuations for T < Tdis- For A > 1 (and 77 = 1/2), 



T, 



dis 



K(T|i3d/eoA2)i/2, whereas for A < 1, Tk 



With the above values this 
= 45^ 0], and Tkis = 83X 



K4(Tdi./eoA)3£oA/ln[7rA/d]. 
results in Tdis = 8/f for Tdi. 
for Tdis = 21QK 0. 

Next we consider the continuous anisotropic case e > 
0. Inspection of the perturbation theory ||] shows, that if 
the disorder induced interaction between FLs is neglected 
the vdW interaction can be written in the following form 



vdW 



x{CUk) + 2CT{k)CdUk)), (7) 
where the interaction between two vortex segments is ^] 



27ri?AVl + A2p 



with Ki, the I'th order modified Bessel function. The 
two differ ent correlations Ct and Cdis are defined as 
Crik) = ( ufcu_fc) - (ufcK^ = T/{ei{k)e) |l| and 
Cdis(fc) = (ufe)(u_fe). The first contribution in (m) re- 



(Ufe)(u_fe) 

suiting from Crik), yields the thermal vdW interaction 
which was considered in Ref. The second contribu- 
tion comes from the disorder and vanishes as the impurity 
concentration vanishes. To find Cdis(fc) we have to dis- 
tinguish the cases Afc <C 1 and Afc 3> 1. For Afc <C 1 we 
may neglect the dispersion of ei{k) and use the known 
result for Ug on large scales ||]. This gives Cdis(fc) ^ 
^2fc"i(A(Afc))2C/3. On the contrary, for Afc > 1, the 
strong dispersion of ei{k) essentially decouples pancake 
vortices in different layers. Hence Cdis ~ Upd, apart 
from logarithmic corrections. Using the asymptotic be- 
havior of the Bessel function Ki{z ^ 0) ~ and 
Ki{z — > 00) ~ e~^, it is easy to see, that the main 
contribution to the second integral in comes from 
large fc, fc ^ tt / {d s R) . In a convenient interpolation 
form the disorder induced van der Waals interaction for 
A < i? < A/e is then given by 

K!S'(«)»-$(Ao^)'(^)'^(.n^)-- 

(8) 



Eqn. (||) is the main result of this paper. For eR <C d 
it changes over to (^). Thus similar to the thermal case, 
the vdW attraction decays as or R~^ for eR < d or 
eR > d, respectively. Eqn. (^) can indeed be written in 
the form of the thermal vdW interaction Q if we replace 
there T by Tdis- 

In the last part we analyze the phase diagram at low T 
and H ^ Hd- The bare repulsion between the flux lines 
2Ko{R/X) and their vdW attraction result in a mini- 
mum of their interaction energy at a distance i?min ~ ctX 
/). For not too low temperatures a is about 20 and 
only weakly T dependent. The same applies to the width 
of the minimum which is of the order of /3A with f3 ~ 10. 
Since the vdW attraction is strongly distance dependent, 
its main contribution comes from those configurations 
where the line pair is at a distance i?min. To lowest order 
in Ao, we can estimate the average vdW interaction by 
considering the configurations of a single line in the ab- 
sence of any FL interaction. With itw A(L/LT,dis)'' for 
the displacement of a single FL, we find from u ~ I for 
the mean distance L|| between two line segments reach- 
ing a minimum i|| « Lx.disCV'^)"'^^^- The length Lg of the 
segment over which the line stays in the minimum follows 
from the same argument as Lg ~ ix.dis P^^'' ■ Thus, the 
contribution of the vdW attraction to the Gibbs free en- 
ergy density is of the order 



7^Kdw(^min)-r^ ~ Kdw(-Rinin) (7 )^'''^^^ 



A2 - 



(9) 



which is much larger than the vdW interaction at the 
mean distance /. For thermal fluctuations (C — 1/2) the 
mean vdW attraction has therefore the same / depen- 
dence as the entropic repulsion. In this case, the result 
can also be obtained by mapping the problem onto 2d- 
Bosons For disorder induced fluctuations, however, 
the vdW attraction decays faster than the steric repul- 
sion. The total Gibbs free energy density for e-Rmin S> d 
can be written in the following form 



G{x;H,T,Tdis) - , 
A 



{zKo(x) -I 5 1- 



dis 



,8/5 



6/5 

hluK} (10) 



which has to be minimized with respect to a: = l/X. 
Here z = 6 is the number of nearest neighbors for 
a triangular lattice. 7t ~ 9.08(r/£oA)2 § and 
7dis ~ Cdis(Tdis/eoA)^'^''^^ denote the strength of the 
entropic and disorder dominated steric repulsion re- 
spectively. Expression (|l^) has to be considered as an 
interpolation between the regimes dominated by the 
bare interaction at high B and the different fluctua- 
tion induced interactions at low B, respectively. The 
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prefactors of the terms following from the vdW inter- 
actions are 6t ~ CT^(l/e)(r/(eoAln^ ^)) and 6dis ~ 
Cdis^(l/e)(rdie/eoA)3''«;-2+6.|i^^|-(i+2.)(^/^)-i+,^ 

The precise determination of the coefficients ct, Cdis, 
Cdis as well as of a and (3 is beyond the scope of this pa- 
per. To find them would require an RG treatment similar 
to that performed in for short range repulsive FLs 
in the thermal case. In this case, instead of singling out 
a typical distance i?min, the contributions of the vdW 
interaction from all distances will be taken into account 
according to their statistical weight. We postpone this 
study to a future publication ||20|| and discuss here the 
phase diagram (Fig. 1) only qualitatively. 

Bragg glass 




fdis Tcu T 



Fig 1. Miikherji and Naltermann 



The coexistence regime is well below the transition to 
the Bragg glass phase where the quasi long range order 
of the lattice persists. This transition line is beyond the 
reach of this analysis and corresponds to a much higher 
value of Bbg 

In conclusion, we have obtained a disorder induced van 
der Waals attraction and a steric repulsion between the 
flux lines in anisotropic or layered superconductors. A 
qualitative analysis of the low field phase diagram for 
such impure system predicts a rich phase diagram with 
a first order transition between two pinned gas/liquid 
phases for not too strong disorder. For larger disorder the 
first order transition disappears and B ^ {H — Hd)^^^ 
decreases continuously as H Hd- 
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FIG. 1. The low- B phase diagram of an impure anisotropic 
superconductor. Thick lines correspond to continuous transi- 
tions. White regions indicate the coexistence regimes of the 
adjacent phases. 5 = corresponds to the Meissner phase. 
For sufficiently weak disorder and low T a low density pinned 
gas phase is separated from a high density pinned liquid phase 
by a first order transition where B jumps from Bi to B2 . With 
increasing disorder, the Bi and B2 lines merge as indicated 
by the dashed line, leading to a critical point at T^. For large 
disorder Td moves towards Tcu and the first order transition 
disappears completely. 

For T = 0, to begin with, we find from ( plj| ) that 
for increasing H the transition from the Meissner phase 
to the pinned liquid phase is continuous since the steric 
repulsion dominates over the vdW interaction for large 
X — l/X and hence B ~ h^^^. Increasing H further, for 
weak enough disorder (but A > 1) the vdW attraction 
will dominate over the steric repulsion in an intermediate 
range of x resulting in a discontinuous transition from a 
dilute to a dense pinned liquid phase. At even higher 
B a second discontinuous transition to the Bragg-glass 
phase takes place. At finite but low T this picture is es- 
sentially unchanged and will then cross over smoothly to 
the thermal phase diagram studied in Q. For stronger 
disorder the steric repulsion dominates at T = for all 
values of B and the discontinuous transition disappears. 
The first order transition between the low (or zero) B 
phase and the dense pinned liquid phase is now shifted 
to higher T and will disappear for sufficiently strong dis- 
order completely. The latter is most likely the situation 
in BSSCO with impurity concentration corresponding to 
Tdis = 4,5K or higher. 
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